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IN  TEODUCTION 

Originally,  the  primary  object  of  thi»  project  wai  the  etudy  of 
contrail  formation  for  altitudes  above  fifty  thousand  feet,  directed 
toward  a  further  knowledge  of  contrail  formation.  More  accurate 
forecasts  could  then  bo  made.  Also  with  a  better  understanding  of 
the  physics  of  contrails  a  method  of  suppression  or  eventual  elimination 
of  contrails  could  be  developed. 

Thd  basic  data  for  this  study  was  to  be  obtained  from  firings  at 
Cape  Canaveral,  Florida.  It  proved  to  be  virtually  imposslbla 
to  obtain  sufficient  data  upon  which  to  base  a  etudy.  While  awaiting 
this  data  it  was  found  that,  although  a  method  of  forecasting 
formation  of  contrails  below  fifty  thousand  feet  had  been  doveloped, 
there  was  still  much  work  to  be  done  in  this  field.  Also,  a  need  did 
exist  for  a  method  of  forecasting  the  height  at  which  contrail*  formed. 
Attention  was  then  focused  on  these  two  problems. 

During  the  period  1  December,  1954  to  15  December,  1955  extensive 
contrail  dhta  was  collected  by  the  U.S.A.  F.  under  Air  Force 
"Project  Cloud  Trail".  This  data  was  made  available  to  Weather 
Services  as  a  basis  for  the  further  study  pf  the  subject  of  contrail 
formation  and  altitude  at  which  these  contrails  develop. 


ABSTRACT 


To  forecast  the  occurrence  ol  contrail*,  a  simple  mathematical 
predicting  model  was  developed.  This  was  based  primarily  oa 
the  shape  of  the  ascent  curve.  The  mathematical  problem  was 
one  of  discriminant  analysis,  in  which  a  linear  combination  of 
three  slopes  was  found,  which  maximized  the  ratio  of  the 
variability  of  the  resulting  combination  between  the  two  categories, 
contrails  and  no  contrails.  The  discriminant  values  found  for  all 
cases  separated  well  ranging  from  85%  to  100%  correct. 

a  4  t 

In  the  forecasting  of  the  height  of  contrhll  formation,  an  equation 
was  developed  using  the  prevailing  height  of  the  tropopauee  to 
indicate  approxiroat  height  of  contrail  formation  and  the  height 
at  which  the  -50<b  t  vnperature  occurred.  The  data  was  normalised 
providing  a  forecasting  method,  effective  regardless  of  location. 

A  graph  of  this  equation  resulted  in  a  simple  method  of  forecasting 
contrail  height. 


DATA  USED 


The  data  used  resulted  from  Air  Force  "Project  Cloud  Trail'*. 

This  study  was  published  in  February  1956  in  an  AWS  Manual 
(AWS  TR  105-132).  The  data  was  collected  from  1,  December, 

1954  to  15,  December,  1955,  in  the  vicinity  o £  twenty-one  upper 
air  sounding  stations.* 

The  procedure  used  in  the  observational  phase  of  "Project  Cloud 
Trail"  was  as  follows: 

a.  Dally  at  approximately  one  hour  before  or  two  hours  after 
1530Z,  two  aircraft  were  vectored  to  point  25,000  feet  above 
an  upper-air  sounding  station.  The  aircraft  then  climbed  to  the 
maximum  altitude  obtainable,  maintaining  position  within  30 
iniles  of  the  station. 

b.  The  wingman  observed  whether  or  not  the  lead  aircraft 
produced  exhaust  trails,  noting  if  they  were  continuous  or 
intermittent,  weak  or  strong,  Including  basis  and  tops  of 
layers  in  which  the  trails  formed. 

/ 

In  this  study  observational  data  for  Caribou,  Maine  and  Sault 

St.  Marie,  Michigan  was  insufficient  to  warrant  their  inclusion. 

.  •-  } 


DAT^.  LIMITATION 

During  the  period  of  data  collection,  the  maximum  altitude  of 
the  aircraft  eeldom  exceeded  45,~000  feet.  Therefore,  this 
study  was  restricted  to  data  observed  below  the  height.  Uelng 
this  forecasting  technique  on  data  for  altitude  above  45,000 
feet,  should  be  accomplished  with  this  restriction  in  mind. 
This  in  no  way  implies  that  the  "altitude  formation  forecast 
method"  developed  herein  would  not  be  valid  at  higher  levels, 
but  merely  that  its  validity  has  not  been  tested. 


FORECASTING  THE  OCCURRENCE  OF  CONTRAIL  FORMATION 

METHOD 

Contrail  data  was  examined  with  the  corresponding  soundings  endeavoring 
to  categorize  the  conditions  necessary  for  formation.  Such  variables  1 
as  tropopause  height  and  lapse  rate  proved  to  be  related  to  the  formation, 1 
therefore  a  method  of  describing  the  general  shape  of  the  temperature 
curve  was  needed. 

In  order  to  completely  describe  any  temperature  ascent  curve.  It 
would  be  necessary  to  use  an  infinite  number  of  points.  However, 
in  this  case,  ouly  that  portion  of  the  ascent  associated  with  the 
existence  of  contrails  required  description.  This  was  accomplished 
by  selecting  values  at  three  set  points.  Generally  speaking,  the 
slopes  of  the  curve  at  -45*C,  -50’C  and  -55*C  were  used  to- 
characterize  the  behavior  of  the  curve  although.  In  order  to  obtain 
the  best  results,  the  three  temperatures  used,  differed  slightly 
from  station  to  station.  These  slopes  were  computed  on  the  basis 
of  the  temperature  change  (AT)  over  a  range  measured  from  one 
thousand  feet  above  and  below  the  three  temperatures  selected. 

(see  table  1)  Thus  three  variables  as  slope  measurements  were  to 
characterize  the  general  behavior  of  the  ascent  curve  In  the  vicinity 
of  the  contrail  occurrence. 

The  mathematical  problem  was  one  of  discriminant  analysis.  In  whidh 
a  linear  combination  of  the  three  slopes  was  found,  so  as  to  maximize 
the  ratio  of  the  variability  of  the  resulting  combination  between  the 
two  categories,  contrails  and  no  contrails,  as  compared  with  the 
variation  within  the  given  category.  A  linear  combination  of  the 
measurements  was  determined  which  attempted  to  place  the  result¬ 
ing  number  in  the  correct  category. 

In  this  discriminant  analysis,  three  variables  were  associated  with 
each  case  of  contrails  and  a  similar  three  variables  with  each  case 
of  no  contrails;  the  three  variables  being  the  slope  of  the  ascent  at 
-45,  -50  and  -55*C.  The  discriminant  process  is  similar  to  a 
regression  analysis  in  whidh  three  are  only  two  categories,  as  In 
this  case  namely  contrails  and  nc  contrails. 
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TABLE  X 


TEMPERATURES  VARIABLES  FOR  PITTSBURG 


CATEGORY  I  (CONTRAILS)  CATEGORY  H  (NO  CONTRAILS) 


-45*C 

-50*C 

-55*C 

-45*C 

-50*C 

-55*C 

-5.6 

,  -5.3 

-5.8 

-4.4 

-2.0 

0 

-4.5 

-4.0 

-3.5 

-3.1 

-1.6 

0 

-4.4 

-4.1 

-3.0 

-2.2 

-2.6 

-1.5 

-5.3 

-4.3 

-3.8 

-4.2 

-4.1 

-2.9 

-6.4 

-5.0 

-4.5 

-3.6 

-0.7 

0 

-5.2 

-4.9 

-4.5 

-3.3 

-1.5 

0 

-6.0 

-5.6 

-4.4 

-3.8 

0 

0 

-5.0 

—4. 9 

-4.0 

-5.0 

0 

0 

-5.9 

-5.0 

-5.4 

-3.5 

-1.8 

0 

-5.2 

-5.0 

-4.5 

-4.0 

-4.2 

-1.6 

-4.8 

-5.0 

-4.5  _ ~ 

-4. 8 

-5. 1 

0 

-5.0 

-5.0 

-5.5 

-4.4 

-1.8 

-1.8 

-4.8 

-5.3 

-4.7 

-2.4 

-1.8 

-3.3 

-5.1 

-5.5 

-5.0 

-3.0 

-3.1 

-5.5 

-5.0 

-4.8 

-0.4 

0 

-4.6 

-5.0 

-5.6 

-1.4 

0 

-5.1 

-5.4 

-5.1 

-3.1 

-0.4 

-2.5 

-5.0 

-4.7 

-4.9 

-4.5 

-5.0 

-1.1 

-4.7 

-4.0 

-3.4 

-4.7 

-5.0 

-5.1 

-5.0 

-4.9 

-4.1 

-6.0 

-5.7 

-6.0 

-3.6 

-4.0 

-3.5 

-4.7 

-5.0 

-4.5 

ANALYSIS 

f 

The  first  step  in  the  analysis  was  to  compute  the  pooled  variance 
and  co-variance  of  these  three  variables  for  each  of  the  categories 
by  adding  the  sum  of  the  souares  and  the  sum  of  cross  products  In 
each  category.  This  was  then  divided  by  the  total  number  of  degrees 
of  freedom,  namely  the  sum  of  the  number  of  observations  In  the 
contrail  category,  plus  the  number  of  observations  In  the  no-contrail 
category  minus  two.  The  mean  of  each  of  the  variables  must  be 
eliminated  so  that  the  sum  of  squares  and  the  sum  of  cross  products 
are  about  a  mean  of  aero. 

The  pooled  variances  and  co-variances  for  the  above  data  are  shown 
in  the  first  three  columns  of  Table  II.  The  fourth  column,  is  the 
difference  between  the  averages  of  variable  -45,-50  and  -55  for 
the  two  situations. 


TABLE  II 


-45  *C 

-50*C 

-55*C 

-45*  C 

0.657675 

v  t 

0.114825 

¥ 

-0.015626 

-1.558 

-50#C 

0.114325 

1.034075 

0.421924 

-2.998 

-55*C 

^■0. 01 5626 

-.421924 

1.322178 

-3.466 

The  values  of  Table  II  are  actually  the  coefficients  of  three  linear 
equations  which  are  solved  for  operators  L/,  L^  and  L^.  The 
variables  are  multiplied  by  these  operators  to  maximize  the  ratio 
of  the  variability  between  the  two  categories,  as  compared  to  the 
variability  within  categories. 


ThMqaaitioB  ffordetwrinfaftht  operators  arc 

(1)  (0.657675)1^+  (0.11 4825) L^+C-0. 015626)1,^*  -1.558 

(2)  (0.114825)1,,  -(1.034975)L^-(0. 4219Z4)L^=-2.998 

(3)  (-0.015626^ -(0.421924)L;t-(1.322178)L}3  i.466 
therefore: 

L, *  2.099912,L^=  1.821953,  L  ?=  2.064852 
substitute  these  values  In  equations  (1),  (2),  and  (3)  above  and 
compute  the  sum  for  each  case. 

The  same  operator  is  used  on  both  the  contrail  and  no-contrail 
data  producing  a  set  of  numbers.  When  good  discrimination  is 
present,  one  category  will  produce  numbers  of  an  order  of 
magnitude  different  from  that  of  any  other.  In  this  study, 
such  a  point  was  computed  which  separated  the  "contrail"  category 
from  the  "no-contrail"  category. 

The  discriminator  values  obtained  in  this  manner  for  Pittsburgh  are 
shown  in  Table  HI.  The  separation  point  was  determined  to  be  twenty 
(20),  on  which  basis  98%  of  the  forecasts  are  correctly  discriminated. 

The  discriminant  values  for  each  station  studied  bad  a  definite 
point  of  separation.  (See  table  IV)  Considering  all  cases,  the 
accuracy  of  the  predictions  of  "Contrail"  vs"No  Contrail" 
ranged  from  85%  to  1C0%.  This  table,  in  addition  to  providing 
the  separation  point  for  each  station  and  the  degree  of  accuracy 
obtained,  gives  the  temperature  points  used  and  the  operators 
computed  for  each  station..  - 

In  utilizing  this  technique  for  forecasting  occurrence  of  contrails, 
three  operators  and  a  separation  point  for  the  discriminant  values 
have  been  determined.  If  "universal  values"  for  both  operators  and 
separation  points  could  he  derived  the  forecasting  technique  would  be 
greatly  simplified. 

To  allow  comparison  of  these  operators,  the  values  were  normalized 
so  that  the  sum  of  the  Ly,  L^,  and  L?  equaled  a  minus  one.  The 
separation  points  were  also  normalized.  The  close  slmiliarity. 
between  stations,  Indicated  that  an  average  set  of  both  operators  and 
discriminant  value  might  very  well  be  developed  which  would  be  valid 
for  all  locations.  This  averaging  yielded  the  following  values: 
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TABLE m 


discriminant  value 

DISCRIMINANT  VALUE 

35.4 

Contrail 

24.0 

Contrail 

33.4 

ii 

22.9 

ii 

32.7 

22.3 

No  Contrail 

31. 8 

•• 

22.1 

Contrail 

31.8 

it 

20.8 

n 

31.1 

•i 

19.4 

No  Contrail 

31.1  - 

*  «• 

19:3 

ii 

31.0 

ii 

17.2 

•i 

10.  & 

•i 

16.2 

ii 

30.3 

•i 

13.1 

ii 

29.5 

ii 

12.9 

n 

29.4 

ii 

12.5 

ii 

29.3 

•• 

12.4 

n 

29.2 

ti 

10.6 

29.1 

ii 

10.6 

ii 

28.5 

ii 

9.7 

ii 

28.3 

ii 

9.4 

ii 

27.9 

ia 

8.8 

ii  * 

27.7 

it 

8.2 

ii 

26.8 

it 

8.1 

ii 

24.2 

ii 

8.0 

ii 
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L,»  0.14 
La*  0.32 
0. 54 

and  a  separation  point  of  2. 8. 

The  effectiveneaa  of  theae  ao  called  "universal  values"  In  fore¬ 
casting  the  occurrence  or  non=occurrence  of  contrails  below  45 
thousand  feet  waa  tested  on  the  data  for  each  station.  It  was  found 
that  the  average  values  were  equally  effective  at  all  stations.  There 
fore,  by  utilizing  this  mathematical  technique  which  describes  the 
slope  of  a  temperature  ascent  curve,  one  can  quite  accurately 
forecast  the  formation  of  contrails. 


FORECASTING  TECHNIQUE 

The  proceedure  to  be  used  is  as  follows:  from  the  latest  raob,  . 

(1)  determlne^T  for  1000  feet  above  and  below  the  levels 
at  which  the  temperatures  -45*C,  -50*C  and  -55*C  occur. 

(2)  From  Table  V,  obtain  the  values  corresponding  to  the 
>  AT  for  L'jL^ind  L ^  (The  purpose  of  thla  table  was  to 

eliminate  unnecessary  multiplication). 

(3)  Total  these  three  values  and  if  this  sum  Is  greater  than 
2.  8,  then  contrails  may  be  expected  to  develop,  but  if  the 
sum  Is  less  than  2.8  then  no  contrails  are  expected. 


os  as  a  Jo  Aj 


TABLE  V- 
(cont) 


L2  x 
-50* 


-T.I 

-7.2 

-7.3 

-7.4 

-7.5 

-7.6 

-7.7  ‘ 

-7.8 


8J 
.2 
.3 
.4 
.5 
.6 
-8.7 
-8.8 
-8.9 


r  »  • 

-9.2 

-9.3 

-9.4 

-9.5 

-9.6 

-9.7 

-9.8 

-9.9 

-10.0 


1.0 

1.0 

1.0 

1.0 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.3 

1.3 

1.3 

1.3 

1.3 

1.3 

1.3 

1.4 
1.4 
1.4 
1.4 


2.3 

2.3 

2.3 

2.4 
2.4 

2.4 

2.5 
2.5 

2.5 

2.6 
2.6 
2.6 
2.7 
2.7 

2.7 

2.8 
2.8 
2.8 
2.8 
2.9 
2.9 
2.9 
3.0 
3.0 
3.0 
3.1 
3.1 

3.1 

3.2 
3.2 


3.8 

3.9 
3.9 
4.0 
4.1 

4.1 

4.2 

4.2 

4.3 

4.3 

4.4 

4.4 

4.5 

4.5 

4.6 

4.6 

4.7 

4.8 

4.8 

4.9 
4.9 
5.0 
5.0 
5.1 

5.1 

5.2 

5.2 

5.3 

5.3 

5.4 


FORE'CASTTI'tC ALTITUDE  OT  CONTRAIT.  FORMA HOIT 


With  a  method  established  of  forecasting  the  existence  or  non 
existence  of  contrails,  the  next  stop  was  to  forecast  the  altitude 
at  which  the  contrails  would  occur.  Indications  were  that  the  base 
of  the  contrail  might  be  directly  related  to  such  meteorological 
variables  (upper  air  parameters)  as  tropopause  height  and 
temperature,  humidity  of  surrounding  air,  turbulence  and  location 
of  jet  stream  core,  lapse  rates  at  various  levels  and  perhaps  the 
relative  location  of  the  general  circulation  patterns  at  constant  pressure 
levels,  say  above  300  mbs. 

•- 

As  mentioned  previously,  "Project  Cloud  Trail"  data  was  used. 
Altitudes  of  all  phenomena  were  reported  in  thousands  of  feet 
according  to  the  ship's  altimeter.  The  reported  altitudes  were 
converted  to  pressure  altitude. 

RAOB's  were  now  plotted  for  all  the  data  cards  for  which  altimeter 
settings  were  available.  Following  discussions  on  properties  affecting 
the  commencement  altitude  of  the  contrail  formation,  it  was  felt  that 
the  one  atmospheric  phenomena  that  most  affected  contrails  was 
the  tropopause.  In  an  attempt  to  determine  how  the  tropopause 
affected  formation,  graphs  of  tropopause  height  versus  contrail 
height,  tropopause  height  versus  contrail  temperature,  tropopause 
temperature  versus  contrail  height  and  tropopause  temperature  versus 
contrail  temperature  were  plotted  for  a  number  of  stations.  Figures 
I  thru  IV  are  the  plots  for  Tuscon,  Arizona  which  are  representative 
of  all  stations  in  this  study. 

Figurel.  The  Temperature  of  the  contrail  base  varies  within  a 
small  area  regardless  of  season;  while  the  tropopause  temperature 
decreases  with  season,  being  colder  in  the  summer  and  warmer  in 
the  winter. 

Figure  n.  The  tropopause  height  varies  with  the  season  while  the 
contrail  temperature  remains  with  in  a  small  range. 

..Tigure  III.  Illustrates  the  relationship  between  contrail  height  and 
tropopause  temperature.  Although  a  relationship  exists,  it  would 
necessitate  taMng  temperatures  from  a  sounding  with  an  unattainable 
degree  of  accuracy.  Therefore,  it  was  felt  that  using  tropopause 
temperature  in  order  to  estimate  contrail  height  was  not  the  answer. 
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Figure-  IV.  The  relationship  between  the  tropopamse  height  end  the 
altitude  at  -which  contrail*  formed  approximated  a  straight  line.  It'* 
slope  was  dependent  ou  the  season  of  the  year.  Figure  IV  became 
the  basis  of  the  first  parameter.  Slmlliar  graphs  were  plotted 
for  all  remaining  Air  Force  Bases  with  the  exception  of  Limestone 
AFB  and  Kinross  AFB  for  which  insufficient  data,  -was  available. 
(Copies  of  all  graphs  are  at  end  of  report). 

Tho  first  step  in  estimating  contrail  height  was  to  locate  the  ' 

prevailing  height  of  the  tropopause  on  Figure  IV,  which  in  turn  will 
indicate  the  approximate  height  of  the  contrail  formation.  This 
proved  quite  sucessful.  However,  no  one  parameter  is  sufficient 
in  itself  to  forecast  heights  of  contrail  formation  therefore  a 
second  parameter  was  needed  which,  when  taken  Jointly  with  the 
first  would  serve  as  a  satisfactory  basis  for  forecasting  the 
altitude  of  contrail  formation. 

An  attempt  was  made  to  relate  the  vertical  location  of  the  jet 
stream  core  to  the  hase  altitude  of  contrails.  It  was  found  that 
vertical  movement  of  the  jet  stream  did  affect  the  contrail  height, 
however  the  tropopause  split's  at  the  jet  core,  wish  the  tropopause 
height  to  the  left  of  the  jet  axis  being  lower  than  cm  the  right. 

In  order  not  to  re-introduce  this  effect  into  any  computations 
the  jet  stream  was  eliminated  from  the  study. 

No  measured  humidity- data  was  available.  However  the  presence’ 
or  lack  of  moisture  has  been  accounted  for  in  the  forecast,  by  the 
beha-dor  of  the  lapse  rate.  A  stable  lapse  rate  inferring  descending 
and  therefore  dry  air  while  an  unstable  lapse  rate  implies  ascending 
and  therefore  moist  air.. 

300znb,  200mb,  and  150mb,  facitnile  maps  were  ct-sed  to  study  - 
the  effects  of  the  laTge  scale  circulation  patterns  -cm  formation. 

This  was  accomplished  by  plotting  the  relative  diecaace  of  tho 
center  of  troughs  in  degrees  longitude  east  or  wresft  of  the  station 
in  question.  It  was  felt  that  such  large  scale  motions  could  not 
be  related  to  the  height  fluctuations  of  contrail  formation. 

A  variable  which  would  reflect  lapse  rate  over  a  period  of  time 
was  worth  investigating.  Contrails,  it  would  appear,  occur  in 
most  cases  at  high  altitudes  and  therefore  at  very  cold  temperatures 


-1*- 


£  revious  studies  have  established  that  wLth  a  temperature  below 

-60  *tT  the  relative  frequency  of  visible  cent  rails  increases  markedly.. 

For  our  purpose  the  height  at  which  -50*C  first  occurs  was  taken  as 

a  separation  point  for  visible  and  non  visible  contrails.  Of  course 

the  bases  of  contrail  formation  do  not  always  occur  at  heights 

above  -50*C,  but  instances  of  bases  occurring  at  warmer  temperatures 

are  generally  explained  by  sudden  drops  in  the  tropopause  height. 

Thus  -50*C  represents  a  limit  below  which  strong  contrails  will 
generally  occur.  The  height  at  -50*C  when  viewed  overs  few  success¬ 
ive  soundings  does  reflect  a  change  In  lapse  rate. 


Our  theory  with  regard  to  the  parameters  affecting  contrail  formation 
was  that  relatively  few  parameters  significantly  effected  the 
formation  while  a  myriad  of  parameters  affect  formation  to  a  minor 
degree.  With  this  in  mind,  further  investigation  of  additional  parameters 
was  halted.  The  two  basic  parameters  were  tested  by  running  a 
linear  correlation. 


(l)Z,-.;+a|Xl+a^  . 

Where: 

Z'  is  the  predicted  contrail  height, 

Xjls  the  estimated  contrail  height  determined  from  Figure  IV. 
X^s  the  height  at  which  -50*C  occurs  bn  the  latest  sounding. 

The  a*  ,  &j  and  a^are  constants  derived  from  the  distribution  of  the 
contrail  data  for  each  air  base. 

In  the  case  of  Pittsburg,  52  cases  of  both  strong  and  moderate 
contrails  were  used  to  determine  the  distribution  of  Z,.  the  contrail 
height  obtained  from  "Project  Cloud  Trail".  The  equations  of  least 
square  are: 

*•“  +  *£x,  +  "Iz 

*£V  »*&iVfzxi 

Their  solution  when  substituted  Into  (1)  yields  the  desired  polynomial 
of  best  fit. 

Table  VI  illustrates,  the  percent  reduction  achelved  for  each  station 
when  this  method  is  tested  using  the  observed  data  from  the  Cloud 
Trail  cards  and  the  RAOBS. 
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TABLE  VI 


„  AIR  FORCE  BASE-UPPER  AIK  STATION  ♦  PER  CENT  REDUCTION  -SAMPLES  SIZE 


KIRKLAND  AFB  (ABQ) 

82.0 

118 

GEIGER  AFB  (GEG) 

94.3 

58 

DAVIS-MONTHAN  AFB  (TUS) 

75.3 

41 

OIXS  AFB  (ACK) 

79.4 

40 

GREAT  FALLS  AFB  (GTF) 

78. 1 

35 

LANGE  LEY  AFB  (NGU) 

73.1 

25 

ELLSWORTH  AFB  (RAP) 

59.8 

94 

McCHORD  AFB  (SEN) 

95.0 

47 

McGUIRE  AFB  (HEM) 

95.8 

50 

MITCHELL  AFB  (HEM) 

94.5 

81 

PORTLAND  NATIONAL  AFB  (PDX) 

91.7 

38 

NIAGARA  AFB  (BUF) 

93.6 

62 

HAMILTON  AFB  (OAK) 

80.5 

44 

ANDREWS  AFB  (DCA)  : 

85.3 

74 

OXNARD  AFB  (LGB) 

86.1 

88 

WRIGHT -PATTERSON  AFB  (FFO) 

84.4 

54 

PITTS3URG  NATIONAL  AFB  (PIT) 

96.4 

52 

GRIFFITHS  AFB  (RME) 

93.6 

44 

SELFRIDGE  AFB  (MTC) 

85.8 

64 

Haring  developed  a  "universal  equation"  for  forecaating  the  oc¬ 
currence  of  contrail*,  it  was  thought  advlskble  to  attempt  a 
■jndllar  “universal,  aquation.11  for  forecaating  the  altitude  at  which, 
contrails  would  develop. 

a»»  and  ivwere  averaged  for  all  stations  and  thir  resulting 
values,  a»  =-ob.  4,  a  =  +  0. 41 ,  and  a^».  +0. 80  were  tested  to 
determine  the  effect  on  the  previously  calculated  percent  reductions 
based  on  constants  derived  for  each  particular  station.  Table  VII 
gives  the  comparison  of  the  original  percent  reductions  to  the  per¬ 
cent  reduction  obtained  from  using  averaged  values.  It  can  be  seen 
that  the  new  method  reduced  the  percent  reduction  to  a  minor 
degree,  namely  about  2%  in  most  cases.  Therefore,  an  averaged- 
equation  can  be  used  to  forecast  the  height  of  contrail  formation 
regardless  of  the  location  of  the  station. 

The  forecasting  proceedure  developed  thus  far  Is; 

1.  From  the  current  sounding  determine  the  altitude  at 
which  -50 *C  occurs.  This  value  is  Xi. 

2.  From  same  sounding  obtain  tropopause  altitude. 

3.  From  figure  IV  obtain  the  estimated  contrail  height 
which  is  then  X|. 

4.  The  altitude  of  contrail  formation  can'then  be  forecasted 
using; 

(1)  Z1  *  a#+  a  |  Xj+'  s^Xj^ 

in  which: 

a,  r-60.4 
-  a,  =  0.41 
a^=  0.  80 

Since  the  relationship  between  the  measured  tropopause  height  (H) 
for  a  particular  day  and  the  estimated  contrail  height  (Xj)  for  that 
day  is  linear,  then 

X/a  e.  +  ®,  H 

Therefore  substituting  this  in  equation  (1)  above  for  each  station 
a,  +  a,e,  +  a^H  +  ^X^ 


results.  Letting 

A  °Ca»  +  aj®.) 
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T3VBLI:  TIT 


Station: 

* 

%  Reduction  Obtained  from 
Value*  for  Individual  Station*: 

.%  Reduction  from  Averaoed 
Value*  for  all  Station*: 

PIT: 

96.4 

95.9 

GEG: 

94.3 

93.3 

NGU: 

73.1 

71.3 

RME: 

* 

93.6  *  • 

92.9 

MTC: 

85.8 

84.5 

OAK: 

80.5 

77.7 

GTF: 

78.1 

76.4 

PDX: 

91.7 

90.8 

SEN: 

95.0 

93.9 

HEM: 

94.5 

93.6 

ADW: 

85.3 

83.3 

WRI: 

96.0 

96.0 

BUFj  . 

93.6 

95,0 

TUS: 

• 

75.3 

73.2 

LGB: 

a*»  -60.4 
'a./  =»  +0. 41 
a^a+0.  80 

86.1 

81.4 

Than:  Z'  »  A+  EH+  CX» 


Averaging  these  computations  for  each  station  resulted  in  a. 
equation  requiring  only  tropopause  height  and  the  height  at  which 
-50*C  occurs  thus  eliminating  the  necessity  of  step  #3  above,  that 
Is,  estimating  contrail  height.  This  average  equation  is: 

Z' »>30.5  +  0.28H  +  0.80X, 

Figure  V  represents  a  graph  of  this  equation  with  the  tropopause 
height  (H)  given  as  the  abscissa  and  the  height  of  the  -50 *C 
(Xz)  given  as  the  ordinate.  The  Intersection  of  H  and  Xg*  Z' 
the  height  in  thousands  of  feet  at  which  contrails  are  lllcely  to  form 
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SUMMARY  AND  EXAMPLE  OF  USAGE 


The  primary  object  of  this  study  was  to  develop  a  mathematical 
method  for  forecasting  both  the  development  of  contrails  and  the 
height  at  which  they  would  form. 

In  forecasting  the  formation  of  contrails,  three  variables  as  slope 
measurement  were  used  to  characterize  the  general  behavior  of  the 
ascent  curve  for  a  particular  location.  Utilizing  these  valuer  dis¬ 
criminant  analysis  was  carried  out  in  which  discriminants  values 
for  each  station  were  determined.  Each  station  studied  had  a 
definite  point  of  separation  between  contrails  and  no  contrails 
discriminants.  The  resulting  operators  and  separation  points  were 
normalized  and  then  averaged  giving  universal  values  for  the 
operators  and  separation  points  which  were  effective,  regardless 
of  location. 

In  forecasting  the  altitude  at  which  contrails  would  occur,  two 
parameters  were  used,  the  height  of  the  -50*C  temperature  and  an 
estimated  contrail  height  based  on  the  pre-established  linear  re¬ 
lationship  between  it  and  the  measured  height  of  the  tropopause. 

The  equation  developed  was: 

JV  »  A  +  BH  +  CXa 
in  which, 

Z1  =  predicted  height  of  the  contrail 
H  =»  Height  of  the  tropopause 
X,  =  height  of  ~50*C 

An  average  equation  for  all  stations  studied  was  determined  to  be: 
Z'  a  -30.5  +  0.20H  +  0.80X, 

This  equation  was  graphed  and  the  intersection  of  H  and  Xs  gave 
the  altitude  of  the  contrail  formation. 

Using  these  two  "universal  equations"  the  following  would  be  the 
proceedura  used  in  forecasting: 

L/  (dT  ©  -45*)  +  L.UT  @  -50*)  +  La(flT  @  -55*  )  *  D 
where  D  =»  separator  point 

0.14  01T  ©  -45*)  +  0.32  (£T  <§  -50*  )  +0.54  ^T  @  -55*)  *  D 
If  D  1b^2.8  no  contrails 

D  le}2.8,  contrails  will  form 


-Z6- 

\  ,  * 


Therefore  from  Table  V 

EXAMPLE:  A.T  at  -45*C  «  -5 

AT  **  “50*C  »  -3 
£T  at  55*C  »  -6 

therefore  the  sum  of  the^T  Is  5. 1  which  la  greater  than  D  so  it 
may  be  concluded  that  contrails  will  form. 

Assuming  from  the  latest  sounding,  thr.t  H  the  tropopause  height 
was  46,000  feet  and  the  X2  the  height  of  the  -50*C  was  40,000  feet, 
by  using  figure  V  we  could  conclude  that  the  contrails  forecasted 
to  form  should  do  so  at  approximately  41, 800  feet. 


RECOMMENDATIONS 


A  simple  but  valuable  forecasting  procedure  has' been  developed 
as  a  result  of  this  study. 

It  Is  therefore  suggested  that: 

1.  this  method  be  tested  throughout  the  country  on 
Independent  data. 

2.  pertinent  data  be  collected  for  altitudes  above  45,000 
feet  to  test  the  validity  of  this  method. 

3.  If,  data  Is  already  available  then, 

a.  this  technique  be  tested  ou  the  new 

b.  changes  be  made  for  Its  Improvement  where  necessary 

c.  efforts  be  made  to  forecast  the  tops  or  endings  of 
the  contrails,  In  this  study  the  endlnga  occurrdd 
above  the  45,000  foot  altitude. 


yivTV-T^> 
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CONCLUSIONS 


If  tiie  forecasting  of  contrails  and  the  commencement  altitudes 
are  still  of  paramount  Importance,  then  the  forecasting  techniques 
developed  as  a  result  of  this  study  should  be  valuable  tools  to 
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ARLINGTON  HALL  STATION;  ARLINGTON  12  VIRGINIA 


NOTICE:  WEE3T  GOVERNMENT  OR  OTHER  EEAWTNGS,  £j?EC  IFICATKL  C-  "  7. 
OTHER  DATA  ARE  USED  FOP.  ANY  PURPOSE  ETHER  THAN  IN  CONNECTION 
WITH  A  DEFINITELY  RELATED  GOVERNMENT PROCURE  jAct  T  i-A  7. 

THEU.  S.  GOVERNMENT  THEREBY  INCURS  NT  RESPONSIBILITY,  RC.Sa2V 
OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE  GOVERNMENT  MAY 
HAVE  FORMULATED,  FURNISHED,  OR  IN  ANTWAY  SUPPLIED  THE  SAID 
DRAWINGS,  SPSEIFICATIONS,  OR  OTHER  DATS.  IS  NOT  TO  B.Z  REGARDED  BY 
IMPLICATION  C2  OTHERWISE  AS  IN  ANY  MANSER  LICE1.  ENG  THE  HOLDER 
OR  ANY  OTHER.  PERSON  OR  CORPORATION,  <CE  CONVEY  NG  ANY  RIGHT3  OR 
PERMISSION  TC  2 AANUFAC TURE,  USE  OR  SELL  ANY  PAT  SN1ED  INVENTION 
THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


